Partonic mean-field effects on matter and antimatter elliptic flows 
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Using a partonic transport model based on the Nambu-Jona-Lasinio (NJL) model, we study the 
effect of scalar and vector mean fields on the elliptic flows of quarks and antiquarks in relativistic 
heavy ion collisions in Au+Au collisions at snn = 7.7 GeV and impact parameter b = 8 fm that 
leads to the production of a baryon-rich matter. Converting quarks and antiquarks at hadronization 
to hadrons via the quark coalescence model, we further study the dependence of the transverse 
momentum integrated relative elliptic flow differences between protons and antiprotons, lambda 
and anti-lambdas, and positively and negatively charged kaons on the strength of the quark vector 
coupling. Our results suggest that a relative weak vector coupling seems to be needed to describe the 
experimental data measured by the STAR Collaboration in the Beam Energy Scan (BES) program 
at the Relativistic Heavy Ion Collider (RHIC). 

PACS numbers: 25.75.Nq, 25.75.Ld 



I. INTRODUCTION 



In relativistic heavy-ion collisions, a hot and dense 
matter that consists of deconfined quarks and gluons is 
produced in the initial stage. Depending on the energy of 
collisions, this so-called quark-gluon plasma (QGP) can 
have various temperatures and baryon chemical poten- 
tials. As the formed QGP expands and cools, it is con- 
verted to the normal hadronic matter (HM). Therefore, 
heavy ion collisions at relativistic energies provide the 
possibility to study the phase structure of the strongly 
interacting matter that is described by the quantum 
chromodynamics (QCD). For the top energy available at 
RHIC and at the Large Hadron Collider (LHC), the pro- 
duced QGP are nearly baryon free and thus has very 
small baryon chemicalpotentials. According to the lat- 
tice QCD calculations the transition from the QGP 
phase to the HM phase in this region of the phase dia- 
gram is a smooth crossover without a clear phase bound- 
ary. The phase transition between QGP and HM is, 
however, expected to change from the crossover to a 
first-order transition at certain finite baryon chemical 
potential called the critical point in the QCD phase di- 
agram (3H(|. To probe this region of the QCD phase 
diagram, the beam energy scan (BES) program at much 
lower energies of v / snn = 7.7, 11.5, and 39 GeV than 
the top energy have recently been carried out at RHIC 
by the STAR Collaboration Q. Although no definitive 
signals for a first-order phase transition and the critical 
end point have been established, a number of interest- 
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ing results have been observed Q. One of them is the 
increasing difference between the elliptic flows of parti- 
cles and antiparticles, thus a breaking of the number of 
constituent quark number scaling of elliptic flows, as the 
collision energy decreases. Such a behavior cannot be 
described by a simple hydrodynamic or hadronic cascade 
model even if the quark coalescence is considered during 
the hadronization of produced QGP Q. Several theo- 
retical attempts have been made to explain this surpris- 
ing experimental result [lol-[l3j|. In particular, the effect 
of hadronic mean-field potentials on the elliptic flows of 
particles and antiparticles has been studied in Ref. [l2[ 
using a multiphase transport (AMPT) model that in- 
cludes both initial partonic and final hadronic scatter- 
ings [Til, [lj|. Because of the different mean-field po- 
tentials for particles and antiparticles in the baryon-rich 
matter formed in these collisions, p, K + and ir~ were 
found to have larger elliptic flows than p, K~ and 7r + , 
respectively, as observed in experiments. However, these 
results are not in quantitative agreement with the exper- 
imental data. The calculated relative integrated elliptic 
flow difference between particles and antiparticles, de- 
fined by [w 2 (particle) — w 2 (antiparticle)]/w 2 (particle), for 
Au+Au collisions at ^/5nn = 7.7 GeV and impact pa- 
rameter b = 8 fm is about a factor of two smaller than 
the measured values of 63 % for p and p and -3 % for 
7r + and 7T _ but is about a factor two larger than the 
measured value of 13 % for K + and K~ . 

The study in Ref. [l2j has neglected the effect of mean- 
field potentials in the initial partonic stage. As shown in 
Ref. [16| using a partonic transport model based on the 
NJL model for Au+Au collisions at ^/snn = 200 GeV, 
where the baryon chemical potential of produced matter 
is small, the attractive scalar mean field on quarks and 
antiquarks was found to reduce their elliptic flows. In 
the present study, we extend the study of Ref. [lj| to 
also include the effect due to the vector mean field in 
the quark matter. Since the latter is repulsive for quarks 
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and attractive for antiquarks in baryon-rich quark mat- 
ter, it increases the elliptic flow of quarks and reduces 
that of antiquarks. Using the coalescence model to con- 
vert quarks and antiquarks to hadrons at hadronization, 
we further show that the vector mean field in the QGP 
also leads to an increase of the elliptic flows of p, A, 
and K~ and a decrease of the elliptic flows of p, A, and 
K + . Our results thus help to explain the large p and p 
but small K + and K~ relative V2 differences observed 
in experiments. Also, information on the strength of the 
quark vector mean field is important for understanding 
the equation of state of QGP at finite baryon chemical 
potential, which at present cannot reliably be obtained 
from lattice studies. As shown in studies based on both 
the NJL and the Polyakov-Nambu-Jona-Lasinio (PNJL) 
models , the existence and the location of the crit- 
ical point in the QCD phase diagram is sensitive to the 
strength of isoscalar vector coupling. A critical point 
only exists if the strength of the vector coupling is small. 
Otherwise, the phase transition is always a crossover. 

This paper is organized as follows: In Sec. UH we 
briefly review the NJL model for three flavors of quarks 
and antiquarks and discuss the mean-field approxima- 
tions. We then describe in Scc. lIIII the partonic transport 
model, with the mean fields taken from the NJL model, 
that is used to study the time evolution of the parton 
phase-space distributions in relativistic heavy-ion colli- 
sions. Results on the quark and antiquark elliptic flows 
from solving the transport model are shown in Sec. IIVI 
together with those of p and p, A and A, and K + and 
K~ that are obtained from quarks and antiquarks by us- 
ing the coalescence model at hadronization. Finally, a 
summary is given in Sec. [V] 



II. THE NAMBU-JONA-LASINIO MODEL 

The NJL Lagrangian for three quark flavors has the 
form M: 



denoting the determinant in flavor space jl8| : 

det/(#Y) = eijkiuTqijidTq^isTqk), (2) 

It gives rise to four-point interactions in two flavors and 
six-point interactions in three flavors. In the two flavor 
case, the sum of scalar and pseudo-scalar interactions and 
the KMT interaction with K=-G reduces to the original 
NJL model @,[l9|. 

Considering only the flavor singlet vector interaction 
in the second interaction term and takingthe mean-field 
approximation, the Lagrangian becomes [H, [2^] 



C = i> (id» - -Gv^i^p^ - i>M*tp + ... (3) 

where M* = diag(M u , M d , M s ) with 

M u = m u - 2G(uu) + 2K(dd)(ss), 
M d = m d -2G(dd) +2K(ss)(uu), 
M s = m s -2G(ss) +2K(uu)(dd), (4) 

and ... denotes constant terms such as (Tpip) 2 . The mean 
fields are calculated as following: 
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7AsE.^m)-fi(k)}, (5) 

where N c is the number of colors, and fi(k) and fi(k) 
are number densities of flavor i and its anti-flavor, re- 
spectively. Because the NJL model is not renormaliz- 
able, the momentum integration requires a cut-off A. For 
the parameters in the NJL model, we use those from 
Refs. @, that is, m u = ra d = 3.6 MeV, m s = 87 
McV, GA 2 = 3.6, if A 2 = 8.9, and A = 750 MeV, unless 
otherwise stated. 
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where ip = (tp u , ip d , tps) T , M = diag(m u , m d , m s ) and A a 
is the Gell-Mann matri ces with A being identity ma- 
trix multiplied by y/2/Z. In the case that the vector and 
axial- vector interactions are generated by the Fierz trans- 
formation of the scalar and pseudo-scalar interactions, 
their coupling strengths are given by Gy = G a = G/2. 
The last term is the Kobayashi-Maskawa-t'Hooft (KMT) 
interaction that breaks U(1)a symmetry with det/ 



III. TRANSPORT EQUATION FOR PARTONIC 
MATTER 



Similar to that for the hadronic matter [23J, |24| , the 
time evolution of the partonic matter produced in rela- 
tivistic heavy ions collisions can be described by the fol- 
lowing transport equation for the partonic phase-space 
distribution function /(x, p); 



If 
dt J 



V x f-V x H-V p f = C, 



(6) 



where ii(x, p) is the hamiltonian of a quark in the self- 
consistent scalar and vector mean fields and C denotes 
the collision term that describes the parton scattering, 
for which we use a constant cross section of 10 mb. 
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For a parton in the scalar and vector mean fields de- 
rived from the NJL model, we have 



H = yJM* 



- p* 



'±gvp = E*±g vP a , 



(J) 



where p* = PT gvP with p = {ipfip) and gy = (2/3)Gy ; ^ 
p° = (ip"f°ip} being the net baryon density, and the upper 
and lower signs are for quarks and antiquarks, respec- 
tively. 

In the test particle method of solving the transport 
equation (25j , their equations of motion are given by 



dxj 
~~dt 
dpi 
dt 



OH _ P * 
dp t ~ E* 

dH 

dx t 

M* dM* 

~ E* dxi 



v dpj _ dpo 
3 dxi dxi 



(8) 



(9) 



We note that in the case that the effective mass of the 
(anti)quark does not change with position, the LHS of 
Eq. © leads to the following Lorentz force: 



Ft 



dt 

= 9v[Vj 



dpi 
dt 

dxi 



Pt 
9v dt 

dp _ dpi 
dxi dt 



= g v (v xB + E)j, 
with B = V x p and E = -Vp° - dp/dt. 

IV. RESULTS 
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FIG. 2: (Color online) Transverse momentum distributions 
of down and strange quarks and their antiquarks in Au+Au 
collisions at ^/snn = 7.7 GeV and impact parameter b — 8 fm 
from the AMPT model. 



For the initial quark and antiquark rapidity and trans- 
verse momentum distributions in a relativistic heavy-ion 
collision, we use the valence quarks and antiquarks con- 
verted from hadrons that are obtained from the Heavy- 
Ion Jet Interaction Generator (HIJING) model [22[ 
through Lund string fragmentation as implemented in 
the AMPT model with string melting [T^|. They are 
given in Fig. [1] and Fig. [2] for down and strange quarks 
and their antiquarks in Au+Au collisions at ^snn = 
7.7 GeV and at impact parameter b = 8 fm. The dis- 
tributions of up and anti-up quarks are slightly smaller 
than those of down and anti-down quarks. It shows that 
the number of down quarks is more than twice that of 
anti-down quarks in heavy ion collisions at this energy, 
indicating that the produced quark matter is baryon-rich. 
We note that strange quarks have a softer transverse mo- 
mentum spectrum than anti-strange quarks because they 
are mostly from the conversion of hyperons while anti- 
strange quarks are mostly from the conversion of strange 
mesons. 

We solve Eqs. (|5|), (JSJ, and © numerically and then 
use the resulting parton momentum distributions to eval- 
uate their elliptic flow vi according to 



''2 



if (Pf 



where N is the total number of test particles. 
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FIG. 1: (Color online) Rapidity distributions of down and 
strange quarks and their antiquarks in Au+Au collisions at 
Y / snn — 7.7 GeV and impact parameter b — 8 fm from the 
AMPT model. 



A. quark and antiquark elliptic flows 

Figure [3] shows the integrated V2 of light and strange 
quarks and antiquarks as functions of time for the cases 
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FIG. 3: (Color online) Integrated elliptic flow «a of light and 
strange quarks and antiquarks at mid-rapidity (\y\ < 1) as 
functions of time in Au+Au collisions at -^/snn = 7.7 GeV and 
at impact parameter b = 8 fm for the cases of including only 
scalar mean field (S), scalar and time component of vector 
mean field (S + Vo), scalar and space component of vector 
mean field (S + Vi), and scalar and both components of vector 
mean field (S + Vo+Vi) 



of including only scalar mean field, scalar and time com- 
ponent of vector mean field, scalar and space component 
of vector mean field, and scalar and both components of 
vector mean field. It is seen that the quark and anti- 
quark V2 are essentially the same if only the scalar mean 
field is included in determining their propagations in the 
quark matter. This is due to the fact that through VM* 
in Eq. © the scalar mean field generates the same at- 
tractive force on quarks and antiquarks. Including the 
time component of vector mean field caused by excess 
quarks leads to a larger quark than antiquark V2 due 
to the resulting repulsive force on quarks and attractive 
force on antiquarks. We note that V2 is generated by the 
gradient of the pressure, which always points outward in 
heavy-ion collisions but is enhanced by a repulsive force 
and reduced by an attractive force. Therefore, the time- 
component of vector mean field enhances the V2 of quarks 
and suppresses that of antiquarks while the attractive 
scalar mean field suppresses both quark and antiquark 
V2 (HI. As to the space component of vector mean field, 
its effect on the V2 of quarks and antiquarks is opposite to 
that from the time component but with a much smaller 
magnitude. The other difference between the effect of 
space component and that of time component of vec- 
tor mean field is that the splitting between V2 of quarks 
and antiquarks starts earlier in the time-component case 
than in the space-component case. The reason is that the 
space component of vector mean field is proportional to 
the baryon current that takes time to develop and is not 
large unless the relativistic flow streams are large [26j |. 
The net result of the partonic mean fields in our study 



is that it leads to a larger quark than antiquark V2 if the 
produced quark matter is baryon-rich. These behaviors 
are seen for both light and strange quarks. On the other 
hand, while the integrated i>2 of strange quarks continues 
to increase with time, that of light quarks decreases at 
later times. This is due to the stronger attractive scalar 
mean field for light quarks than for strange quarks as a 
result of the larger decrease of the light quark conden- 
sate along the radial direction than that of the strange 
condensate and the smaller light quark mass than the 
strange quark mass. 
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FIG. 4: (Color online) Elliptic flow vi of light and strange 
quarks and antiquarks at mid-rapidity (\y\ < 1) as functions 
of transverse momentum at hadronization for the cases of in- 
cluding only scalar mean field (S), scalar and time component 
vector mean-field (S + Vo), scalar and space component vector 
mean field (S + Vi), and scalar and both components of vector 
mean-field (S + Vo + Vi) 

For the transverse momentum dependence of quark 
and antiquark V2 , they are shown in Fig. [4] for those at 
the end of the partonic phase, which is about 2.1 fm/c 
after the start of the partonic evolution when the energy 
density in the center of produced quark matter decreases 
to 0.8 GeV/fm 3 , again for the cases of including only 
scalar mean field, scalar and time component of vector 
mean field, scalar and space component of vector mean 
field, and scalar and both components of vector mean 
field. For the last case, the difference between the inte- 
grated V2 of up and down quarks and their antiquarks is 
about 49 % of the integrated «2 of up and down quarks, 
and that between strange quarks and anti-strange quarks 
is about 34 % of the integrated v 2 of strange quarks. 



B. hadron elliptic flows 

To study how different quark and antiquark V2 is re- 
flected in the i>2 of produced hadrons, we use the coa- 
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lesccncc model to convert them to hadrons at hadroniza- 
tion. In this model, the probability for a quark and an 
antiquark to form a meson is proportional to the quark 
Wigner function of the meson with the proportional con- 
stant given by the statistical factor gM for colored spin- 
1/ 2 quark and antiquark to form a colorless meson (271 — 
[29I1. that is 



/ M (p,k p ) = 8g M exp 



-£--kV 



(12) 



factors in the quark coalescence are = 9k- =1/36 
and g p = g\ = gp = g\ = 1/108. The oscillator fre- 
quency uj is determined from the root-mean-square ra- 
dius of the produced hadron, which is taken to be 0.6 
fm for K+ and K~ , and 0.877 fm [13 for both p and A 
as well as for p and A. We do not consider tt + and ir~ 
elliptic flows because they are affected similarly by the 
partonic mean fields as a result of similar elliptic flows 
for u and d quarks as well for u and d antiquarks. 



where 



1 , \ | rx m 2 k! - mik 2 

P=-=(ri-r 2 ), k p = V2 . 

V2 mi+TO 2 



(13) 



with rrij, and being the mass, position and mo- 
mentum of quark (antiquark) i, respectively. The width 
parameter a p in the Wigner function is related to the 
root-mcan-square radius of the meson via 



2 (mi 
3 



l tomim2(mi + m 2 ) ' 



(14) 



where the second line follows if we use the relation a „ = 
l/y/juicj in terms of the oscillator frequency w and the 
reduced mass [i\ = 2(l/mi + l/m 2 ) . 

The probability for three quarks or antiquarks to co- 
alescence to a baryon or an anti-baryon is similarly pro- 
portional to the quark Wigner function of the baryon, 



f B (p, A,k p ,k A ) 



8~g B exp 



k 2 a 2 



k 2 rr 2 
K X cr X 



(15) 



where gs is the statistical factor for three colored spin- 
1/2 quarks to form a colorless baryon, and 
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m 2 r 2 



- r 3 



kx = 



mi + m 2 
3 m 3 (ki + k 2 ) - (mi + m 2 )k 3 
mi + m 2 + m 3 



(16) 



The width parameter a\ is related to the oscillator fre- 
quency by (/z 2 o;) _1 , with [i 2 = (3/2) [1/ (mi + m 2 ) + 
l/m 3 ] _1 . The root-mean-square radius of a baryon or 
an antibaryon is then given by 
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FIG. 5: (Color online) Elliptic flows V2 of midrapidity (\y\ < 
1) p and p (left panel), A and A (middle panel), and K + and 
K~ (right panel) at hadronization as functions of transverse 
momentum. 

Figure [5] shows the w 2 of p and p (left panel), A 
and A (middle panel), and K + and K~ (right panel), 
at hadronization as functions of transverse momentum. 
It is seen that the quark coalescence leads to a larger 
hadron u 2 than the quark i> 2 at same transverse mo- 
mentum. Furthermore, the v 2 of p, A, and K~ are re- 
spectively larger than those of p, A, and K + , leading 
to the relative differences between their integrated u 2 , 
[v2 (particle)— v 2 (antiparticle)]/v 2 (particle), of 36.5, 44.3, 
and -5.0 % respectively, compared with the values 63±14, 
54±27, and 13±2 % measured in experiments. 

The dependence of the relative difference between inte- 
grated particle and antiparticle v 2 on the vector coupling 
gv is shown in Fig. [S] Besides the value gv/G = 0.33, 
which is from the Fierz transformation, two other values 
of 0.165 and 0.73, the latter from Rcf. [2l| . are also used. 
The relative integrated v 2 differences are seen to increase 
almost linearly with the strength of the vector coupling 
but become saturated beyond the value gv/G = 0.33. 



1 m\(m 2 + m 3 ) + m 2 (m,3 + mi) + m|(mi + m 2 ) 
2 



w(mi + m 2 + m 3 )mim 2 m 3 



(17) 



For the mesons K + and K , the baryons p and A, and 
the anti-baryons p and A considered here, their statistical 



V. SUMMARY 

We have studied the effect of partonic mean fields on 
the elliptic flows of quarks and antiquarks in a brayon- 
rich quark matter by using a transport model based on 
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FIG. 6: (Color online) Relative differences between integrated 
V2 of mid-rapidity (\y\ < 1) p and p, A and A, and K~ and 
K + at hadronization for several values of the isoscalar vector 
coupling. 

the NJL model. For the scalar mean field, which is at- 
tractive for both quarks and antiquarks, it leads to a 
similar reduction of the quark and antiquark V2 as first 
found in Ref. [Hj]. The vector mean field, on the other 
hand, has very different effects on quarks and antiquarks 
in the baryon-rich matter as it is repulsive for quarks and 
attractive for antiquarks. The time component of vector 
mean field turns out to have the strongest effect, result- 
ing in a significant splitting of the quark and antiquark 
i)2 as a result of enhanced quark vi and suppressed anti- 
quark V2- The space component of vector mean field has, 
however, an opposite effect; it suppresses V2 of quarks and 
enhances that of antiquark, although relatively small and 
appearing later in the partonic stage compared to that 
of the time component of vector mean field. 

Using the quark coalescence model, we have further 
studied the elliptic flows of p, A, and K + and their 
antiparticlcs produced from the baryon-rich quark mat- 
ter and found that the differences between particle and 



antiparticlc elliptic flows are appreciable as a result of 
the different quark and antiquark V2- The magnitude 
of the relative integrated V2 difference between particles 
and their antiparticles depends on the strength of the 
vector coupling. For the one obtained from the Fierz 
transformation of the scalar interaction in the NJL model 
(g 1) /G=0.33), the V2 differences between p and p, A and 
A, and K + and K~, are about 36.5, 44.3, and -5.0 %, 
respectively. Compared to the data from the BES ex- 
periments at RHIC, the calculated value for the relative 
V2 difference is about a factor of two smaller for p and 
p, slightly smaller for A and A, and about a factor of 
two smaller but with an opposite sign for K + and K~ . 
Since hadronic mean fields can further increase the rel- 
ative V2 differences between p and p, A and A, and K + 
and K~ [l2T |. a smaller vector coupling in the baryon- 
rich quark matter seems to be needed. Also, allowing 
chemical reactions of partons such as pair creation and 
annihilation of quark and antiquark is also expected to 
lead to a splitting of the quark and antiquark V2 as shown 
in a recent study based on the ideal hydro+UrQMD hy- 
brid model [l| , which assumes local thermal and baryon 
chemical equilibrium in the initial stage of heavy ion col- 
lisions. A more quantitative study that includes both 
the partonic and hadronic mean-field effects as well as 
the chemical reactions in the partonic phase is needed 
before a definitive determination of the strength of vec- 
tor interactions in the baryon-rich QGP can be achieved. 
Nevertheless, the upper bound inferred from the present 
study for the value of the vector coupling already has 
an important implication that a critical point is likely 
to exist in the QCD phase diagram based on the NJL 
model. 
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